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m . 1 INTRODUCTION 

o- 

rvj ■ How matter is distributed around galaxies is one of the 

r- H fundamental questions in cosmology. Gravitational lens- 
ing provides a powerful method to map the matter dis- 
trib ution at small and large lengt h scales (for a review, 

S^ see iBartelmann fe Schneider! 120011 ). Recent large galaxy 
•_h . redshift sur veys, such as t he Sloan Digital Sky Sur- 

_C^_. vey (SPSS; lYork et all l200Ch and COSMOS survey, have 
allowed us to explore the mean surface density pro- 
file of galaxies through weak lensing techniques (e.g., 
Sheldon et alj|2004 IMandelbaum et alfaod : IMenard et all 
2010 j; ILeauthaud et al.ll2012T i. 

IMenard et"a l. (2010a, hereafter MSFR) measured the 
mean surface matter density profile of the SDSS main galax- 
ies with the mean redshift (z) = 0.36 through gravitational 
lensing magnification of background quasars (QSOs). They 
calculated the cross-correlation between the number density 
of foreground galaxies and the flux magnification of back- 
ground QSOs. The cross-correlation function was then con- 
verted to the surface matter density profile £ m of the lens 
galaxies as a function of the projected distance R from the 
galactic center. The derived profile is well approximated as 
E m oc R-° 8 at lOkpc < R < lOMpc. 

MSFR also detected the systematic offset between the 



ABSTRACT 

We develop an analytic halo model for the distribution of dust around galaxies. The 
model results are compared with the observed surface dust density profile measured 
through red dening of background quasars in the Sloan Digital Sky Survey (SDSS) 
reported by IMenard et al.l ( 2010al) . We assume that the dust distribution around a 
galaxy is described by a simple power law, similarly to the mass distribution, but 
with a sharp cut-off at aR v [ r where i? v ir is the galaxy's virial radius and a is a model 
parameter. Our model reproduces the observed dust distribution profile very well 
over a wide range of radial distance of 10 — 10 4 /i _1 kpc. For the characteristic galaxy 
halo mass of 2 x 10 12 /i _1 Mq estimated for the SDSS galaxies, the best fit model is 
obtained if a is greater than unity, which suggests that dust is distributed to over a 
few hundred kilo-parsccs from the galaxies. The observed large-scale dust distribution 
profile is reproduced if we assume the total amount of dust is equal to that estimated 
from the integrated stellar evolution over the cosmic time. 
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five SDSS photometric bands in magnification of background 
QSOs; in shorter wavelength, QSOs appear less magnified. 
It is interpreted as reddening due to dust in and around 
foreground galaxies. Adopting the small Magellanic cloud 
type dust model for the sample galaxies, they derived the 
mean surface dust density profile of galaxies T,d(R) from the 
galaxy-QSO color cross-correlation function. The shape of 
the derived £d is very similar to that of E m at lOkpc < R < 
lOMpc, suggesting that there is a substantial amount of dust 
in the galactic halos (see also Chelouche et al. 2007; McGee 
& Balogh 2010). Theoretical models are needed to properly 
interpret the observationally inferred dust distribution. 

In this Letter, we develop an analytic model based on 
the so-called halo approach to study the distribution of dust 
aroun d galaxies. Earlier in iMasaki. Fukugita fc Yoshidal 
(|2012l ). we used cosmological A-body simulations to study 
in detail the matter distribution around galaxies. There, 
we showed that the observed surface density profile can be 
used to determine the characteristic mass of the sample lens 
galaxies, and that the mass distributed beyond the galax- 
ies' virial radii contributes about half of the global mass 
density. We provide a physical model for the dust distri- 
bution in this Letter. Our model is characterized by two 
key physical parameters; one is the host halo mass of the 
galaxies and the other is the extent of dust distribution. 
The for mer is determined from the observed matter pro- 
file (e.g.. IMandelbaum et al.l l200rj; lHavashi fc White! 120081 ; 
ILeauthaud et al.ll2012l ; IMasaki. Fukugita fc Yoshidal 12012? ). 
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while the latter can be inferred from the dust distribution. 
How far dust is transported from the galaxies is indeed a 
highly interesting question. The observed dust profile is well 
described by a single power law over a wide range of distance 
of from lOkpc to lOMpc. We show that the profile is decom- 
posed into two parts, the so-called one-halo and two-halo 
terms. We parametrize the one-halo term such that dust is 
distributed to aR v [ r where R v [ r is the galaxy's virial radius. 
Through model fitting, we determine the extension param- 
eter a to be greater than unity. We discuss the implication 
for the dust production and transport mechanism into in- 
tergalactic space. 



2 THE MODEL 

2.1 Halo approach 

We present a simple formulation to calculate the surface dust 
density pr ofile. Our model is based on the s o-called halo 
approach jSeliakl l200(il ; ICoorav fc Shethl 120021 ). The mean 
surface density Ed(-R) is divided into two terms: 

Z d (R) = T.f(R) + Zf(R), (1) 

where R is the distance in the projected two-dimensional 
plane. The one-halo term T, d h (R) arises from the central 
halo, and the two-halo term E d h (i?) from the neighbouring 
halos. 

The contribution from an individual galaxy halo with 
mass Mh to the one-halo term E d h (7?) is given by the pro- 
jection of the halo dust density profile pd(r\Mh) along the 
line-of-sight \: 



Y, d {R\M h ) = 



dxPd(r=^ X 2 + R 2 \M h ). 



(2) 



The one-halo term is then a number-weighted average of 

X d (R\M h ) 
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(4) 



where An/ AMh is the halo mass function and M m i n is the 
threshold halo mass for the sample galaxies. The threshold 
mass corresponds to the typical host halo mass of the ob- 
served galaxies. 

We calculate the two-halo term power spectrum Pj h (fc) 
as follows: 



PT{k) = 



Plm(fe) 

_L_ 

Pd 

1 

Mhalo 
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An 
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M d {M h )b{M h )u d {k\M h ) 



AM h 
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(5) 



where p d is the mean cosmic dust density, Pu n (k) is the lin- 
ear matter power spectrum, b(Mh) is the halo bias factor, 
M d (Mh) is dust mass in and around a halo with mass Mh, 
and u d (k\Mh) is the Fourier transform of the density pro- 
file p d normalized by its dust mass. The power spectrum is 
converted to the two-point correlation function via 



*™ = 2^' dkk 



> sin(fcr) 
kr 



P d {k). 



(6) 



Then we obtain the two-halo term of the mean surface den- 
sity profile 



Vf{R) 



= Pd 



2p d 



d X fd h (r=y / X 2 + R 2 ) 



Ar 



rtf 



Vr 



R 2 



(-) 



We adopt a flat-ACDM cosmology, with fi m = 
0.272, »a = 0.728, H = 70.2km s" 1 Mpc _1 ,n s = 0.961 
i|Komatsu et alj201lh . We use t he code CAMB to obtain the 
linea r matter powerspectrum (|Lewis. Challinor fc Lasenbvl 
20001) and utilize the ha lo mass function and bias given by 
Sheth fc Tormenl (|1999l ) at z = 0.36 which is equal to the 
mean redshift of the galaxy sample used in MSFR. 



2.2 Dust distribution profile 

We assume that the spatial distribution of dust within and 
around a halo is organized as 



Pd(r\M h ) oc 



■ exp 



V C*-Rvir / 



(8) 



where R v i v is the virial radius. Within the virial radius 
i?vir, the mea n internal matter d ensity is A x p cr i t , where 
A is given by iBrvan fc Normanl (|1998T l . Essentially, we as- 
sume that the dust distribution follows a singular isothermal 
sphere (SIS) profile with exponential cut-off at r = aR vlT . 
One of our aims in this Letter is to determine the value of a, 
i.e., how far dust is distributed from galaxies. Fig. [T] shows 
the shape of the dust density profile pd for a halo with mass 
Mh = 10 13 /i _1 Mq computed in the above manner. We see 
the dependence on a clearly. 

The power-law shape is motivated by the fact that the 
observationally derived surface dust profile itself is well fit by 
a simple power law of Ed oc R~°' s , similarly to matter distri- 
bution (MSFR). Also, detailed ca lculations of dust e j ection 
and radiation-driven transport bv lBianchi fc Ferraral (|2005f ) 
show approximately a power-law distribution for the result- 
ing gas metallicity through dust sputtering. We have also 
examined other profiles of the form r~ 3 and r -1 with a sim- 
ilar exponential cut-off. However, we have found that nei- 
ther of the steeper or the shallower profile reproduces well 
the observed dust profile at small distances. We therefore 
adopt the profile equation ((8| in our model. 

The Fourier transform of p d (r) is given by 



u d {k\M h ) = 



Ar 4nr' 



sin(fcr) p d {r\M h ) 
kr M d (M h ) 



(9) 



Note that the value of u d should be unity in small-fc limit. 
We determine the amplitude of p d by setting the halo dust 
mass associated with a halo to be a certain value Md- To 
this end, we first consider th e total amount of dust around 
galaxies in the local universe. iFukugital (|2011f ) estimated the 
total amount, in units of the cosmic density parameter, 



^galaxy dust = 4.7 X 10 



(10) 



Interestingly, this value is close to the difference between the 
estimated amount of dust produced and shed by stars over 
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Figure 1. The model dust density profile as a function of the 
spatial distance from the center. The black, red and blue lines 
represent the model profiles with a = 0.1, 1 and 10, respectively . 
For comparison, we also show the NFW IINavarro et al,l [1997]) 
and the untruncated SIS profiles by grey lines. Note that the 
amplitudes are arbitrary in this figure. 



the age of the universe an d the summed amount fo und in lo- 
cal galactic discs (see also llnoue fc Kamavall200of ). Suppose 
that the comoving density of the total halo dust remains 
constant in the local universe. Then the mean cosmic den- 
sity of dust in galactic halos is given by 

pd(z) = figalaxy dust Pclit(z)(l + z) 3 . (11) 

We set the dust mass associated with a halo to be 



dr4nr 2 p d (r\M h )=M d = FxM h (12) 



where F is the dust-halo mass ratio. We integrate the dust 
mass weighted by the halo mass function to obtain the global 
dust density. We normalize pd, or equivalently V, by match- 
ing the global dust density to equation (|ll[l . Note that F 
is not necessarily a constant but can be a function of halo 



2.2.1 Dust-halo mass ratio F 

An essential physical quantity in our model is the dust-halo 
mass ratio F in equation (|12[) . We propose two simple mod- 
els. The first one is constant model, i.e., F is independent 
of halo mass. The dust-halo mass ratio is simply the global 
density ratio 



1.73 x 10" 



n 



galaxy dust/^m 



/fin 



(13) 



Because the heavy elements that constitute dust are 
produced by stars, it may be reasonable to expect that the 
dust mass is proporti onal to the stellar mass. Intriguingly, 
iTakeuchi et~al1 (|201Ch used data of AKARI and GALEX 
to show a moderate correlation between the stellar mass 
and dust attenuation indicator for the sample galaxies (see 
their Figure 16). In our second model, we consider the ob- 
served galaxy stellar-halo mass relation to model the halo 
mass dependence of F. We call the model as mass dependent 



Figure 2. Two models for dust-halo mass ratio as a function of 
halo mass. The dashed and the solid lines are the ratio V for the 
constant and the mass dependent model, respectively. 



model. iLeauthaud et al.l (|2012f ) recently studied the stellar- 
halo mass relation from the joint analysis of galaxy-galaxy 
weak lensing, galaxy clustering and galaxy number densities 
using the COSMOS survey data. We use their functional 
form with the best fit parameters atz« 0.37, 



log 10 (M fc ) = 



log 10 (Mi)+/31og 10 
(Af„/Af., )' 



+ 



l + (M,/M,,o) 



M. 



0.5, 



(14) 



where M* is the galaxy stellar mass, log 10 (Mi/Mo) = 
12.52, log in (M», / M B ) = 10.92, /3 = 0.46,£ = 0-57 , and 
7 = 1.5 (see also iBehroozi. Conrov. fc WechsleiJl2010T ). We 
then relate the dust mass to the stellar mass as 



M d (M h ) oc M*(M h ) 



(15) 



The normalization constant is determined to be 3.05 x 10 
by integrating this equation weighted by the halo mass 
function. The global dust mass density thus calculated is 
matched to equation (|lip . 

Fig. [2] compares V for our two models. The shape of F 
for the mass-dependent model reflects the stellar-halo mass 
relation. The peak value of T at ~ 6 x W 11 h^ 1 Mq is ~ 
10 -4 . Overall, F for the mass dependent model is larger 
than that for the constant model at the characteristic mass 
of the sample galaxies (see Section 3). 

We are now able to calculate the dust surface density 
profile. In Fig. [3] we show the dependence of the surface dust 
density profile on the extension parameter a. For this figure, 
the threshold halo mass is fixed to be 2 x W 12 Ii~ 1 Mq. We 
compare three cases; a = 0.1, 1 and 10. The dotted lines 
show the one-halo term. Clearly the extension parameter a 
affects significantly the amplitude and the shape of the one- 
halo term. The central surface density is larger for smaller 
a. This can be easily understood by noting the total dust 
mass associated with a halo is given by equation (|12|l . On 
the other hand, a does not affect much the two-halo term at 
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Figure 3. The surface dust density profile as a function of the 
projected radius for a = 0.1, 1 and 10. The projected radius is the 
physical distance at z = 0.36. The constant model for dust-halo 
mass ratio is adopted. The dotted and the dashed lines represent 
the one-halo and the two-halo terms, respectively. The solid linos 
show the sum of the two terms. 



Figure 4. The mean surface dust density profile as a function of 
the physical projected distance from the center of galaxies at z = 
0.36. The results from our constant model and mass dependent 
model arc shown by the green and the black lines, respectively. 
The dotted, the dashed and the solid lines are the one-halo, the 
two-halo and the total, respectively. 



R > IMpc. The amplitude of the two-halo term is essentially 
set by the halo bias 6(M n ) (see equation [5]). 



3 RESULTS 

We fit our dust distribution model to the observed surface 
dust density through least chi-squared minimization. We 
have two physical parameters, M m i n and a, in our model. 
We have found that poor constraints are obtained on the pa- 
rameters when both of them are treated as free parameters. 
Because M m j n is already estimated to be 2 x 1O 12 /i _1 M0 in 
iMasaki. Fukugita fc Yoshidal 1)20121 ) through detailed com- 
parison of the observed matter distribution with the results 
of TV-body simulations, it is sensible to fit our dust distribu- 
tion model by treating only a as a free parameter. Namely, 
the characteristic halo mass can be determined by the grav- 
itational lensing measurement of the matter distribution, 
whereas the dust distribution extension can be inferred from 
the observed dust profile. 

We evaluate the likelihood of the specific model by the 
X 2 value of the model fit to the observed quantities. The 
obtained best-fit a for the constant and the mass dependent 
models are, respectively, 



-, 1C +0.203 /-, \ 
= 1-16-0.155 (Iff) 



= 2.. 



-,+0.450 
'-0.355 



(Iff) 



for constant model, (16) 

for mass dependent model. (17) 



Fig. U shows the best-fit dust profile of the constant model 
with a — 1.16 and that of the mass dependent model 
with a = 2.88. The data points are from MSFR. Both 
models for F reproduce the observed profile fairly well. It 
is interesting to compare these two equally good models. 
The mass dependent model requires a larger a, which is 
owing to the difference in the typical value of V for the 
two models. At M.% > 10 1 h~ Mq, F of mass dependent 
model is higher than that of constant model. Because the 



one-halo term is largely contributed by halos with masses 
~ M m i n = 2 x 10 12 /i _1 M Q , the best fit a is larger for the 
mass dependent model to match the observed inner dust 
surface density profile (see Fig[3]). 

Overall, our simple models reproduce the observed dust 
profile very well. Intriguingly, both our models suggest 
a ~ 0(1), i.e., halo dust is distributed over a few hun- 
dred kilo parsecs from the galaxies. It is also important 
that the observed power law surface density Ed oc R~ ofi at 
R = lOkpc - lOMpc can be explained with a ~ 0(1). The 
apparent large-scale dust distribution is explained by the 
two-halo contributions. Dust is distributed to/over ~ 7? v ; r 
from a galaxy, but not necessarily up to 10 Mpc as one might 
naively expect from the observed dust profile. 

It is worth discussing the total dust budget in the uni- 
verse. The amplitude of the two-halo term depends largely 
on the mean cosmic density of intergalactic dust, pd in equa- 
tion q The excellent agreement at large separation be- 
tween the observed dust density profile and our model pre- 
diction shown in Fig. [4] implies that f2 ga iaxy dust should be 
~ 10 -6 . Clearly, a significant amount of dust exists around 
the galaxies. Such "halo dust" can close the cosmic dust b ud- 
get as discussed more quantitatively bv lFukugital (|201lr ). 

The intergalactic dust could cause non-negligible ex- 
tinction and thus could compromise cosmological stud- 
ies wit h distant supernovae (jMenard. Kilbinger fc Scranton I 
uOlObT ) . We calcul ate the mean ext inction by the intergalac- 
tic dust following IZu et al.l (|2011f ) (see their equation [2]). 
With our model mean cosmic density of f2 ga iaxy dust = 
4.7 x 10~ 6 , the predicted mean extinction is (Ay) — 
0.0090 mag to z — 0.5. Such an opacity is not completely 



The halo bias 6(M n ) is also a critical factor. However, 
the characteristic halo mass, and hence fe(M n ), is well con- 
strained from the observed matter density profile, as shown in 
IMasaki. Fukueita fc Yoshidal J2012ft 
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negligible even in the current SNe surveys, and will be- 
come important for future surveys that are aimed at deter- 
mining cosmological parameters with s ub-percent precision 
JMenard, Kilbinger fc Scranton ll2010bT ). 
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4 DISCUSSION 

We have shown that our halo model can reproduce the 
dust profile around galaxies measured by MSFR. By fit- 
ting the model to the observed dust profile, we infer 
that dust is distributed beyond the virial radius of a 
galaxy. Several authors proposed radiation-driven transport 
of dust from galactic discs into intergalactic medium a t 
high reds h ifts JAguirre et alj|200ll ; iBianchi fc Ferrarall2005l ) . 
IZu et al.l (|201lr ) showed that galactic winds can disperse 
dust into the inter-galactic medium efficiently. Such stud- 
ies generally suggest that dust can travel up to a few x 100 
kpc from galaxies if the ejection velocity is ~ 100km s . 
The relatively larger extent radius of dust for our mass de- 
pendent model requires very efficient transport mechanisms. 
Note also that the dust must survive on its way through the 
galactic halos. Dust in a large, group-size halo could be de- 
stroyed by thermal sputtering in hot gas (Bianchi & Ferrara 
2005; McGee & Balogh 2010). Clearly, detailed theoretical 
studies on dust transport are needed. 

Fluctuations of the cosmic infrared background (CIB) 
provide insight into dust distribution a r ound galaxie s (e.g. , 
IViero et al.ll2009l ; lAmblard et "all 1201 ll ). IViero et al.l (|2009h 
used BLAST data to measure the CIB power spectrum. Us- 
ing a halo approach, they found that the observed power 
spectrum at small angular scales is reproduced if halo dust 
extends up to a few times of the virial radius of galactic ha- 
los. It is remarkably consiste nt with our conclusion in this 
Letter. lAmblard et al.l (|201lh compared their measurements 
of t he CIB anis o tropie s from Herschel wide-area surveys 
with IViero et al.l (J2009I ). Two power spectra are consistent 
with each other at small scales. Our dust distribution model 
may provide a key element for studies on the CIB. 

Although our model reproduces the observed dust pro- 
files very well, a few improvements can be certainly made. 
The dust extension a and the dust-halo mass ratio T are 
likely to depend on the halo mass and galaxy type etc 
(McGee & Balogh 2010) . One may need to consider the 
distribution of satellite galaxies within a halo by using, 
for example, the halo occupation distribution (HOD). In- 
deed, we see slight discrepancies between the model predic- 
tions and the observation in the dust profiles at ~ IMpc 
(Fig. [4|, where the contribution from satellite galaxies 
are non-negligible (fo r more detailed modeling, see e.g., 
iLeauthaud et al.l I2012I ). In principle, the HOD parameters 
can be inferred from the lensing magnification measurement 
presented by MSFR. However, in order to derive the pa- 
rameters accurately, one needs to use additional information 
from observations of th e galaxy-galaxy correlation function 
(JLeauthaud et al.ll2012l ). Including these effect in our model 
is beyond the scope of this Letter, but will be needed in order 
to interpret data from future wide-field galaxy surveys. 
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